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Abstract 

In a recent paper Graham et al constructed oscillating and static universe models which are 
stable with respect to all classical perturbations. Here we show that such universes are quantum- 
mechanically unstable and can collapse by quantum tunneling to zero radius. We also present 
instantons describing nucleation of oscillating and static universes from nothing. 



I. INTRODUCTION 



It has been recently shown [T] that the spacetime of an inflationary universe is necessar- 
ily past-incomplete, even though inflation may be eternal to the future. All past-directed 
timelike and null geodesies, except maybe a set of measure zero, reach the boundary of the 
inflating region of spacetime in a finite proper time (finite affine length, in the null case). 
This indicates that inflation must have had some sort of a beginning. One possibility is that 
the universe could have spontaneously nucleated out of nothing [2]. 

Unlike earlier singularity theorems, the theorem of Ref. fP^ does not rely on Einstein's 
equations and does not assume any energy conditions. To show the incompleteness of a given 
geodesic, all it requires is that the expansion rate averaged along the geodesic is greater than 
zeroQ 

Hav > 0. (1) 

This is a rather weak condition, but it points to a possible loophole in the argument: the 
universe could be static in the asymptotic past. Then Hav = and the theorem does not 
apply. For example, a closed Friedmann-Robertson- Walker spacetime 

rfs2 _ ^^2 _ q2^^^^^2 ^2) 

with a scale factor 

a(t) = ao(l + e^"*) (3) 

and ao,Ho = const describes an inflating universe at t > and is geodesically complete both 
to the future and to the past. The idea that the universe could have started as a static, 
closed space in the asymptotic past has been widely discussed in recent years, under the 
name of "emergent universe" scenario (see, e.g., [3H5] and references therein). 

Construction of emergent universe models is a challenging task. First of all, it is not 
easy to arrange for a static universe to be stable. It is well known that Einstein's static 
model, describing a closed universe filled with nonrelativistic matter and positive vacuum 
energy, is unstable with respect to small perturbations of the radius. Even if the universe is 
initially perfectly fine-tuned, it will be destabilized by quantum fluctuations and will either 
start inflating or collapse to a singularity. Such a universe cannot survive for an infinite 

^ The expansion rate H is defined in terms of a comoving congruence; see [T] for details. 
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time. However, radial stability can be achieved in models with modified gravity [3H5] or 
with "exotic" matter [HHH] (see also U^)- 

A particularly simple model of the latter kind was recently discussed by Graham et al 
[8] . It describes a closed universe with a negative cosmological constant and a matter source 
with equation of state 

P = wp. (4) 
The energy density of the universe is then 

p = A + Poa"='(^+"') (5) 

with A < and po > 0? and the Friedmann evolution equation is 

•2,1 SttG 2 /«\ 
a + 1 = pa . (6) 

The model is radially stable, provided that w satisfies 

-!<«;< -1/3. (7) 

For a perfect fiuid source, the speed of sound Cg can be found from = dP/dp = w. With 
the equation of state this gives < 0, indicating instability with respect to short- 
wavelength compressional perturbations. Hence, it is important that the exotic matter 
source should not be a perfect fiuid [Ej|^ It could, for example, be an assembly of randomly 
oriented domain walls, in which case [11] w = —2/3 and > 0. For this choice of w, Eq. ([s]) 
takes the form 

p = A + Poa-\ (8) 
and the evolution equation ^ has a simple oscillatory solution 

a = a;"^(7 - V7^^cos(wt)) (9) 

where 



^=\/fG|A| (10) 



and 



^ A perfect fluid source could be acceptable if one allows an equation of state P{p) more general than 
Eq. (lil). All one needs is that w = P/p satisfies Eq. (Fzl) and dP/dp > 0. 
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Graham et al in Ref. [8] focused primarily on this special case, which they referred to as 
"simple harmonic universe". 

A static universe solution is obtained from ^ by setting 7 = 1; then a = l/cu. It has 
been shown in [8J that this solution is stable with respect to arbitrarily small perturbations, 
including all scalar and tensor modes. 

Apart from stability, emergent universe models need a mechanism that would trigger 
inflationary expansion after an infinitely long stationary phase. A possible mechanism has 
been suggested in [3]. It involves a massless scalar field with a self-interaction potential 
V{(f)), such that V — )■ const at — — oo. In the stationary regime, the field "rolls" from 
— oo at a constant speed, <p = const. Inflation is triggered when arrives at the non-flat 
region of the potential. A similar mechanism has been employed in Ref. 

At the classical level, the simple harmonic universe model, supplemented with a suitable 
mechanism to trigger inflation, yields a consistent emergent universe scenario. Our goal in 
this paper is to investigate whether or not the model remains stable in the quantum theory. 
We shall restrict the analysis to the simplest minisuperspace model with a single dynamical 
degree of freedom - the radius of the universe a. We shall see that already at this simplest 
level the emergent universe exhibits a quantum instability. 

II. COLLAPSE THROUGH TUNNELING 

We consider a spherical universe ^ with a matter content described by Eq. ([s]) . The scale 
factor a (if:) is the single dynamical degree of freedom. In classical theory, such a universe 
can be regarded as a constrained dynamical system with a Hamiltonian 



n 



G 



(12) 



37ra 



where 




aa 



(13) 



is the momentum conjugate to a and 




(14) 



The Hamiltonian constraint 7/ = then yields the evolution equation ml 
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In quantum theory, the universe is described by a wave function ipi^.): the conjugate 
momentum pa becomes the differential operator —id/ da and the constraint is replaced by 
the Wheeler-DeWitt (WDW) equation [12] (for a review see, e.g., [T3HT5] ) 

•H^ = 0, (15) 



or 



rf2 



(3 d 



+ U{a) i){a) = 0. 



(16) 



da'^ a da 

Here, the parameter /3 represents the ambiguity in the ordering of the non-commuting fac 



tors a and pa in the Hamiltonian (12). Its value does not affect the wave function in the 



semiclassical regime a ^ Ipianck- From here on we set /3 = 0. 

One might expect that for a simple harmonic universe the potential U{a) should be of 
the same form as for a harmonic oscillator. This, however, is not the case: the motion in 



the potential (14) is simple harmonic only for a particular value of the energy, = 0. With 



p(a) from (pj), we have 



U{a) 



i) -'^ 



{poa + Aa^ 



(17) 



It will be convenient to introduce a rescaled variable x = coa with u) from Eq. (10). In 
terms of this variable the WDW equation takes the form 

d^ 



dx^ 



+ U{x) ^(x) = 0, 



where 



U{x) = A"V(1 -27x + x^), 



7 is given by Eq. (11), and 



A 



16G^|A| 
9 



(18) 



(19) 



(20) 



The classically allowed range is defined by U{x) < 0. This range is non-empty when 
7 > 1. The shape of the potential in this case is illustrated in Fig. 1. In the classical 
solution, the radius of the universe oscillates forever between the values x+ and x_ where 
f/(x±) = 0, 

x± = 7 ± V72 - 1. (21) 



However, it is clear from the figure that quantum-mechanically the universe can tunnel 
through the barrier to a vanishing size at x = 0. The WKB tunneling action is given by 



S= ^U{x)dx (22) 
Jo 

and the corresponding tunneling probability can be estimated as 

V ~ (23) 

This can be interpreted as the probability of collapse through quantum tunneling as the 
universe bounces at radius x = x_. 

Semiclassical quantum tunneling in oscillating universe models has been studied by 
Dabrowski and Larsen p5]. They considered a closed universe containing nonrelativistic 
matter (dust), a domain wall fluid with equation of state w = —2/3, and a negative cos- 
mo logical constant. Due to the presence of dust, this model has another classically allowed 



range at small values of a. The WKB action (22) for tunneling between the two classi- 
cally allowed regimes can then be expressed in terms of elliptic integrals. In the absence of 
dust, the model of p2] reduces to the simple harmonic universe, but the authors have not 
discussed this case. 



For a simple harmonic universe, the integral in (22) can be expressed in terms of elemen- 
tary functions. 



Since the tunneling probability (23) is nonzero, such a universe cannot survive forever. 

For 7 = 1, the classically allowed range reduces to a single point, and the WKB action is 
given by the simple formula 

S.... = ly (25) 



The classical solution in this case is a static universe with a; = 1, and Eq. (23) can be 



interpreted as being proportional to the probability of quantum collapse per unit time. 



III. TUNNELING FROM NOTHING 



We note that the tunneling between x = X- and x = can also go in the opposite 
direction, in which case Eq. (23) with S from (24) or (25) can be interpreted as describing 



U(x) 




FIG. 1. The WDW potential U{x) for the parameter values A = 0.03 and 7 = 1.1. 

spontaneous creation of an oscillating or static universe from nothing. The corresponding 
instanton can be found by solving the Euclideanized Friedmann equation, 

= uj'^{x+ — x){x^ — x), (26) 

where the dot stands for differentiation with respect to the Euclidean time r. The solution 

can be expressed as 

UT = r = -21n f v^^Z^±V^^ _ (27) 



'0 ^/{x+-x'){x^ -x') V V^+V^ 

Solving this for function of r we find 

x{t) = 7 - ^(7 - - lil + l)e--^ (28) 



Introducing 

m = uj~^ In I 

.7-1 



ro = u-Hn(^], (29) 



Eq. (28) can be rewritten as 



x{t) = 7 — v/7'- lcosh[a;(r-ro/2)]. (30) 
Note that this is related to the Lorentzian solution (|9]) by a simple analytic continuation, 



as one might expect. The instanton solution (30) starts with x = at r = 0, grows until 



it reaches a maximum value x{tq/2) = x_, and then returns to x = at r = tq. It is 
symmetric with respect to the point t = tq/2. 
The geometry of the instanton, 

ds"^ = dr^ + a\T)dQl, (31) 

is similar to a 4-dimensional ellipsoid. We note that 

a(0) = -a(ro) = 1, (32) 

which indicates the absence of conical singularities. In other words, the "poles" at r = 0, tq 
are rounded off. 



For 7 = 1 the instanton solution (28) simplifies to 



x{t) = 1 - e-'^^ (33) 

It interpolates between a; = at r = and a; = 1 at r — t- oo. The geometry of this instanton 
is that of a cigar. It is rounded off at a = and asymptotically approaches a static sphere 
at large r. The instanton action in this case is given by 

\Sinst\ = ^ y dra [d^ + (coa - if] = —. (34) 



Of course it is the same as in Eq. (25 ). Note that the action is finite, even though the instan- 
ton has an infinite 4-volume. We note also that the boundary term, which is proportional 
to the normal derivative of the boundary volume, vanishes for this instanton. 

Even though there are no conical singularities, a closer examination shows that somewhat 



milder singularities are still present at the poles^J The scalar curvature for the metric (31) 
is 

R = 6a~^(l — — ad). (35) 

The first two terms in the parentheses cancel out at the poles, but in the last term d(0) = 
—'yu 7^ 0, and thus R oc oc r~^. This singularity is integrable, so the instanton action is 
finite. 

It is possible that the curvature singularity can be removed by modifying Einstein's 
equations or the equation of state at small values of a. We could imagine, for example, that 



We are grateful to Jaume Garriga for pointing this out to us. 



8 



for a gas of domain walls the equation of state parameter gradually changes from w = —2/3 
to w = —1 as we approach a = (so the equation of state becomes that of the symmetric 
vacuum in the wall interiors). This would cure the singularity. 

The situation here is somewhat similar to that with the Hawking- Turok (HT) instanton 
|17j . which was proposed to describe quantum creation of open universes. Garriga has shown 
that this singular instanton can be regulated with a suitable matter source [18] and can also 
be obtained by dimensional reduction from a regular instanton in a higher-dimensional 
theory (TH]. It is possible that the Euclidean solutions presented here can similarly be 
regarded as approximations to instantons of a more fundamental theory. 

An important difference between HT and our instantons is that in the HT case the vicinity 
of the singular point makes a significant contribution to the action. For our instantons the 
contributions of singular points are negligible. This indicates that the instanton action and 
the tunneling probability are not sensitive to short-distance modifications of the theory. 



IV. THE WAVE FUNCTION 



Having studied the semiclassical tunneling of the universe, we shall now examine solu- 



tions of the WDW equation (18) for the wave function of the universe ijj{a). By analogy 
with a quantum harmonic oscillator, one might expect the wave function to oscillate in the 
classically allowed range and to decay exponentially in the two classically forbidden ranges 
on both sides of it. However, the situation we have here is rather different. In the case of 
an oscillator, we solve the Schrodinger equation 

^-^+^V^^(x) = i?^(x) (36) 

with boundary conditions iplx — ?■ ±oo) = 0. Solutions exist only for certain values of the 
energy, E = {n + ^)u); this determines the energy spectrum of the oscillator. 

Now, in our case the eigenvalue of the WDW operator is fixed: it is equal to zero. If we 
impose boundary conditions requiring, e.g., that ip^x ^ oo) = ip{x = 0) = 0, the system 
would be overdetermined and no solutions would exist, except for some special values of the 
parameters A and 7. For generic values of the parameters, we have the freedom to impose 
only a single boundary condition. A natural choice appears to be 

ij{x ^ 00) = 0. (37) 



U(x),.A(x) 




FIG. 2. Solution of the WDW equation for the parameter values A = .05 and 7 = 1.3 (red dashed 
line). The WDW potential is also shown (blue line). 



This fully specifies the solution. In the classically forbidden region < x < x_, the wave 
function is a superposition of exponentially growing and exponentially decreasing solutions. 
The solution that grows towards a = will dominate, unless the parameters of the model 
are fine tuned to suppress its contribution. Some numerical solutions to the WDW equation 



(18) are illustrated in Figs. 2 and 3. 



The interpretation of these solutions is not completely clear, since we do not have a well 
established procedure for extracting probabilities from the wave function of the universe (see, 
e.g., Ref. [T3j and references therein). But a nonzero value of V'(O) signals a non- vanishing 
probability of collapse and appears to be inconsistent with the picture of an eternal oscillating 
or static universe. 

Here, we assume that hitting the singularity at x = is fatal for the universe. It is 
conceivable that wave functions similar to those in Fig. 3 could describe an eternal universe 
tunneling back and forth between a finite radius a = uj~^ and a Planck-size nugget. However, 
analysis of this possibility would require a full theory of quantum gravity and is beyond 
our present level of understanding. Our simple minisuperspace model certainly becomes 
inadequate at a ~ Ipianck- 

10 
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FIG. 3. Solution of the WDW equation with A = 0.1 and 7 = 1 (red dashed line). The WDW 
potential is also shown (blue line). The classical solution in this case is a static universe at x = 1. 



V. DISCUSSION 



Our analysis in this paper indicates that oscillating and static models of the universe, even 
though they may be perturbatively stable, are generically unstable with respect to quantum 
collapse. Here we focused on the simple harmonic universe with matter content described 
by Eq. ([s]), but we expect our conclusions to apply to a wider class of models. In particular, 
one could investigate the quantum stability of braneworld, loop quantum cosmology, and 
other modified gravity inspired models]^ 

Is it possible to save the simple harmonic universe from quantum collapse? One possibility 
is to impose the boundary condition 

i;{0) = 0. (38) 
(This boundary condition was introduced in [12]; for a recent discussion see [21j.) Together 



with the boundary condition at infinity (37), this will enforce a relation between the pa- 



rameters of the model 7 and A. As Figs. 2 and 3 illustrate, the value of ip{0) can be either 
positive or negative. This is determined by whether iIj{x) is growing or decreasing near 



Some relevant discussion of quantum cosmology in Horava-Lifshitz gravity models can be found in Rcf. 
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X = x_, which is in turn determined by the number of oscillations N of ijj that fit into the 
classically allowed range a;_ < x < x^. Suppose for definiteness that we decrease A while 
keeping 7 fixed. This makes the potential well deeper, so monotonically increases and 
oscillates between positive and negative values, making one oscillation as changes 
by AA^ ~ 1. By continuity, tp{0) should go through zero twice per such oscillation. Values 
of A ^ 1 correspond to the semiclassical regime, where N ^ 1 and the boundary condition 



(38) can be satisfied by a relatively small change in A. 

Thus, for each value of 7 > 1 we expect an infinite set of values of A for which the condition 



(38) can be enforced. Fig. 4 shows the wave function for a universe with the parameters fine- 
tuned in this way. This approach appears to avoid the collapse, but the following argument 
indicates that it may not be possible to extend it beyond minisuperspace. 



The WDW equation ( 15 ) can be interpreted as stating that the energy of a closed universe 
is equal to zero. Quantum states with different occupation numbers of matter particles have 
different energy of matter, but this energy is exactly compensated by the negative energy 
of gravity, so the total energy is zero. Then one expects that transitions between different 
states should be possible, as long as they have the same conserved quantum numbers. For 
example, there seems to be nothing to prevent spontaneous nucleation of particle-antiparticle 
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pairs. This seems to suggest that the universe will evolve to a state with large occupation 
numbers and high entropy. In terms of the wave function, we expect to be a superposition 
of states with different occupation numbers. The value of ipi^) cannot be fine-tuned for all 
of them. Hence, we expect that quantum collapse cannot be prevented by fine-tuning in 
more realistic models including a quantum matter field. It would be interesting to study 
this issue quantitatively. 

A somewhat puzzling aspect of the WDW equation is that the wave function of the 
universe is independent of time. Following DeWitt [I2], we can interpret this as indicating 
that time should be identified with some semiclassical variable characterizing the universe. 
In other words, clocks, being a part of the universe, should also be described by the wave 
function of the universe. Our minisuperspace model has a single dynamical variable a. If we 
use a to measure time, the model has no other variables whose evolution we can describe as 
a function of time. Moreover, the radius of the universe a is a rather poor clock in static or 
oscillating models. The role of a clock in such a universe can be played by the homogeneous 
mode of a massless, minimally coupled scalar field, The field equation for (p in the 

metric ^ yields 

a^(p = const. (39) 

This shows that classically (j) changes monotonically with time, which makes it a good clock. 
More generally, a clock can be defined whenever the model has some semiclassical variables 
described by a WKB factor in the wave function, as discussed, e.g., in fi2[ [T5t I22ti2^ . With 
such variables included, it may be possible to quantitatively define such concepts as the 
tunneling probability per oscillation period or per unit time. We hope to return to some of 
these issues in subsequent work. 
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